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Abstract 

We describe a simple construction of Kingman's coalescent in terms 
of a Brownian excursion. This construction is closely related to, and 
sheds some new light on, earlier work by Aldous [3] and Warren |23j . 
Our approach also yields some new results: for instance, we obtain the 
full multifractal spectrum of Kingman's coalescent. This complements 
earlier work on Beta-coalescents by the authors and Schweinsberg [7]. 
Surprisingly, the thick part of the spectrum is not obtained by taking 
the limit as a — ► 2 in the result for Beta-coalescents mentioned above. 
Other analogies and differences between the case of Beta-coalescents 
and Kingman's coalescent are discussed. 
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1 Introduction and results 



Let (IT, t > 0) be Kingman's coalescent. This is a Markov process taking its values 
in the set V of partitions of N = {1, . . .}, such that initially 11(0) is the trivial 
partition composed exclusively of singletons, and such that each pair of blocks 
merges at rate 1. Kingman's coalescent was introduced in 1982 by Kingman in his 
seminal paper [15]. A fascinating mathematical object in its own right, Kingman's 
coalescent is also a cornerstone of mathematical population genetics, serving as 
the basic model for genealogies in a well-mixed, selectively neutral population. 

In this paper we provide a new, simple construction of Kingman's coalescent 
in terms of a Brownian excursion. We apply this construction to the study of some 
fine properties of this process, paying particular attention to its behaviour near 
the beginning of time, when microscopic clusters coalesce into a finite number of 
macroscopic ones. This phenomenon is known as coming down from infinity, and 
we are able to describe the precise multifractal spectrum corresponding to clusters 
of atypical sizes at small times. Our construction is closely related to (and, in 
some sense, almost implicit in) earlier work by Aldous [3], Warren [23] and others 
(see section [2] below). However, there are important differences, which will also be 
discussed in Section [2] (for instance, the aforementioned application could not have 
been deduced from these works). This paper complements the picture developed 
in [7] and [5] on the relation between coalescents, continuum random trees, and 
Fleming- Viot type particle systems. 

We now describe our construction of Kingman's coalescent. Let (Bt,t > 0) be 
an excursion of Brownian motion conditioned to hit level 1. That is, let v denote 
Ito's excursion measure on the space of continuous excursions fi* := U£>o^C> 
where 

O c := {/ : [0, C] -> K continuous, f(x) = x £ {0, C}}- 

Let then (B tl < t < Q be a realization of sup s>0 B s > 1). (With a slight 
abuse of notation, we may consider B to be a function defined on [0, oo) by declar- 
ing B(s) = for all s > £.) Let {L(t, x)}t>o, x >o denote a jointly continuous version 
of the local-time process of B, and define 

Z x := L((,x), x>0. 

Thus Z x is the total local time accumulated at level x by the excursion (B s , s > 0). 
Define a P-valued process (II u ,0 < u < 1) as follows. Consider the set {ej}^ 1 of 
excursions of B above level 1, ordered according to their height: that is, 

supei(s) > sup £2(5) > • • • 

s>0 s>0 

Now, fix < u < 1, and consider the set of excursions {ek}i<k<N above level u 
reaching level 1, where iV = N(u) is the number of such excursions. Note that for 
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Figure 1: Construction in Theorem 1. In this picture, we have: Tl u = 
({1,2}, {3, 4, 5,6}). 



every i > 1, e% belongs to exactly one excursion et for some 1 < k < N, and let 
<f> u (i) = k G {1, . . . , N} be this index. Then define II U by declaring that for every 
i,j > 1; i and j are in the same block of U u if and only if <fi u (i) = (j> u (j), that is, 
if and only if £j and Ej fall into the same excursion for some k < N. Our main 
result states that II U is in fact a time-change of Kingman's coalescent. 

Theorem 1. The process (Hjj^,t > 0) has the same law as Kingman's coalescent, 
where for all t > 0, 



1. As the reader has surely guessed, the ordering of the excursions by 
their height is not crucial to this result, but is in the spirit of our use of the 
Donnelly-Kurtz lookdown approach (see section [3] and [7]). 

2. The time-change {U(t),t > 0} satisfies the following properties: U(Q) = 1, 
lim^oo U(t) = 0, and U is continuous and monotone decreasing. 




(1) 



Remark 2. 
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As promised at the beginning of this introduction, we now give some applica- 
tions of Theorem Q] to the study of the small-time behaviour of Kingman's coales- 
cent. Let (n t ,t > 0) be Kingman's coalescent, and let 

F(t) = lim —#{1 < i < n : i is in the same block as 1 in Kt}- (2) 

n— >oo n 

F(t) is the frequency of the block containing 1 at time t, and the existence of the 
almost sure limit in ([2]) (for all times simultaneously) follows from general theory 
on coalescent processes and exchangeable partitions: see, e.g., Proposition 2.8 in 
[8]. In the same manner, one may define Fi(t) for all i > 1 to be the asymptotic 
frequency of the block containing i at time t, thus F(t) = Fi(t) in ([2]). The first 
corollary gives us the behaviour of the typical block size near time zero. While 
this result is well known, our proof is new. Along the way we also provide an 
alternative path to a result of Aldous concerning the asymptotic of the number of 
blocks with a given size. 



Corollary 3. As t — > 0, 



t 



where — > stands for convergence in distribution, and where E and E' are two 
independent exponential variables with parameter 1. 

Our second application is, to the best of our knowledge, new. It concerns the 
existence and the number of blocks (in the sense of Hausdorff dimensions to be 
specified below) with atypical sizes as t — > 0, that is, blocks of size of order i 7 with 
7^1. It turns out that for 7 < 1 (i.e., for anomalously large blocks) we have 
to look at a more precise, logarithmic, scale as there are almost surely no blocks 
whose size will be i 7 for small t. In particular, the sizes of the smallest and largest 
blocks at small time t > are identified. 



For this result, the framework introduced by Evans [13] is very convenient. 
Consider a random metric space (S,d), defined as follows. Define a distance d on 
N by declaring that for every i,j > 1, d(i,j) is the time at which the integers i and 
j coalesce, and let S be the the completion of N with respect to d. It can be shown 
that to every ball of radius t > say, corresponds a unique block of the partition 
Kt. The space S is thus naturally endowed with a measure 77 such that for every 
x £ S and for every t > 0, r](x,t) := r](B(x,t)) is the asymptotic frequency of the 
block of Kt associated with B(x,t). In this setting, the question mentioned above 
becomes: are there points x £ S such that r](B(x, t)) is approximately t 7 as t — > 0, 
and if so, what is their Hausdorff dimension? Define for 7 > 1 

C I \ J r- a V lQg(??(x,t)) \ 

Sthin(7) = <x e S : hmsup — = 7 } . (3) 
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This set corresponds to points of S with atypically small rj(x,t). For thick points, 
consider for all 8 > 0, 

SthickW = {xeS: hmsup ^4 = 4 • ( 4 ) 
[ t^o t\logt\ J 

Theorem 4. 

1. If0<8<l then 

dimS , thick (6>) = 1-6, a.s. 

When 6 = 1, Schick (0) 7^ almost surely, but if 8 > 1 i/ien Schick ($) is empty 
almost surely. 

2. If 1 < 7 < 2 i/ien 

7/7 > 2 tfien SWn(7) 

It may be deduced from the above result that the Hausdorff dimension of S is 
equal to 1 almost surely, a result which was first proved by Evans [13] . This result 
should be compared to Theorem 5 in [7] which describes the multifractal spectrum 
for Beta(2 — a, a) coalescent, a E (1,2). Kingman's coalescent is a limit case for 
this family and corresponds formally to the case a = 2, since the Beta(2 — a, a) 
distribution converges weakly to a Dirac mass at when a — > 2. Interestingly, 
only the "thin points" side of the spectrum is obtained by taking a limit as a — > 2 
in that result: at the power-law scale, the "thick points" part of the spectrum is 
empty, although the limit as a — > 2 exists and is non-degenerate. For instance, 
here the smallest block turns out to be of order t 2 as t — ► 0, while the largest block 
is of order £log(l/i). In the case of Beta-coalescents with parameter 1 < a < 2, 
these quantities are respectively t a /( a_1 ) and t l / a , which does not coincide with 
tlog(l/t) when a — > 2~ . 

We emphasize that the proof of Theorem 2] is in some sense purely conceptual: 
all the work consists of identifying Evans' metric space with a certain Galton- 
Watson tree (which here turns out to be the Yule tree), equipped with its branching 
measure. Theorem |4] then follows automatically from the results of Peter Morters 
and Nan-Rueh Shieh |18U 19j on the multifractal spectrum of these measures. These 
results themselves rely on elegant percolative methods introduced by R. Lyons [T7] . 

2 Previous results and discussion 

We review some earlier results concerning the relation between Kingman's coales- 
cent and Brownian processes. We start by discussing the ideas contained in [4] 



dimS t hm(7) = 1, a.s. 

7 

= a.s. but Sthin(2) 7^ almost surely. 
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and [23] , which are most directly related to the representation of Theorem [TJ 

2.1 Excursions conditioned on their local time 

The notion of continuum random tree (CRT) was developed by D. Aldous in his 
seminal papers [HE], in which a careful treatment of the correspondence between 
excursions and continuum random trees is given. Particular attention is paid to the 
normalized Brownian excursion and to the tree that it encodes through this corre- 
spondence. Early effort bore on the identification of the law of the tree spanned by 
a finite number of leaves chosen suitably at random. Given an excursion function 
e E fi* and p G N, let t±, t%, . . . , t p G (0, C) be pairwise distinct times, and define a 
planar labelled tree with edge-lengths T = T p (e, t\, . . . , t p ) as follows. 

- T contains a root at height and p leaves, with the height of the k th leaf 
being e(t k ). 

- The path from the root to the leaves j and k (1 < j ^ k < p) splits at height 
inf{e(s),s G (tj,t k )}. 

- At each branch point, branches are labelled "left" and "right". 

It is further assumed that all branch points are distinct (i.e., branch points are 
binary). Thus T has p — 1 branch points. Le Gall [16] (see his Theorem 3) first 
identified the distribution of T p (e,t±, . . . ,t p ) on the space of trees when e is a 
normalized Brownian excursion and the t k are independent and uniform random 
variables on [0,£(e)] (and further gives the conditional distribution of e condi- 
tionally given the tree T). Later on, Aldous [3] gave a decomposition result for 
this distribution by identifying the conditional distribution of T given the local 
time profile (£ x ,x > 0) where £ x is the total local time accumulated at x by e. 
This conditional distribution is constructed from a certain inhomogeneous (i.e., 
time-dependent) coalescent process whose main ingredient is Kingman's coales- 
cent. More precisely, this process is defined as follows. First, the height of the 
leaves (xi, . . . , x p ) are i.i.d. samples from £ x dx (which is a probability distribution 
because £(e) = 1 almost surely). Then, thinking of time as running from x = oo 
down to x = and the heights Xi as being the birth times of particles, the law of 
T is such that as x decreases, the particles that are present merge pairwise, each 
independently at rate A/£ x . In other words, the conditional distribution of T given 
I can be thought of as an inhomogeneous coalescent with immigration rate i x and 
coalescing rate 4/£ x . This result may be regarded as the continuum random tree 
counterpart to Perkins' disintegration theorem in the context of superprocesses 
|21j . The proof of Aldous is based on a discrete approximation to the continuum 
random tree, but soon afterwards Warren [23] gave two alternative, direct proofs 
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of Aldous's result. The tools used in those arguments are closely related to earlier 
work by Le Gall [16] and Warren and Yor [24] . 

It thus comes as no surprise that one can embed Kingman's coalescent into 
a Brownian excursion. However, we emphasize that our construction is rather 
different in that the tree which we consider is spanned by vertices at distance 1 from 
the root, rather than by leaves selected uniformly at random. Moreover, it seems 
difficult to use the description of T given above to deduce results about Kingman's 
coalescent. On the other hand, as the reader will see, here these conclusions will 
follow in a straightforward fashion once Theorem Q] is proved. Finally, we believe 
that the computations leading up to Theorem [1] are new, and hope that the reader 
will find them interesting in themselves. 

2.2 Analogies and differences with Beta-coalescents 

The present paper complements earlier results of [6] [7J, which focuses on small-time 
properties of Beta -coalescents with parameter a G (1,2). Beta-coalescents are a 
family of A-coalescent (i.e., coalescents with multiple collisions), where the measure 
A is the density of a Beta(2 — a, a) random variable. (The interested reader may 
consult [7J and references therein.) Kingman's coalescent constitutes a formal limit 
case for the Beta-coalescents as a — > 2 since the Beta(2 — a, a) converges weakly 
to the Dirac mass in zero as a — > 2. It was proved in [7J that for all 1 < a < 2, 
Beta-coalescents can be embedded in continuum stable random trees associated 
with a-stable branching processes rather than Brownian motion. (An excellent 
introduction to continuous random trees can be found in |10j). When a = 2 the te- 
stable branching process is the Feller diffusion, which is closely related to Brownian 
motion. It is therefore natural to suspect a relation between Kingman's coalescent 
and Brownian motion. 

However, we emphasize that there is an essential difference between Theorem 
[1] here and Theorem 1 in [7] . The analogue of Theorem 1 in [7] is the following. 
Let (B s , s > 0) be a reflecting Brownian motion and let t\ = inf{t > : L(t, 0) > 
1}, where L(t,x) is the joint local time process of B. Let v > be such that 
v < sup s<T1 B(s), and for all < u < v, define a partition in exactly the same 
way as in the construction given above Theorem [IJ except that where we used the 
level 1 we now use level v: thus, here the excursions (£i)£i which we consider are 
those above level v (instead of 1) and the e& are the excursions above u that reach 
v. Therefore, using this notation, the partition U u defined for Theorem [TJ is simply 

(with, however, the difference that here B is not a single excursion). 

Now, for all t > 0, let 



V(t) :=mf\s>0: / AZ~ l du > t 
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where Z x = L{t u x) for all x > 0. (Note that V(t) < sup s < n B(s)). Fix T > 
and for < t < T, consider the partition 

fr — u v(T) 
lit ■— L1 v(T-t)- 

We can prove the following result: 

Proposition 5. (n s ,0 < s < T) has the same law as (k s ,0 < s < T), Kingman's 
coalescent run for time T. 

Observe that the main difference between Proposition and Theorem 1 is that 
here the tree is spanned by vertices at a random distance V(T) from the root of 
the continuum random tree, while these vertices are at a deterministic distance 
1 in Theorem [TJ This was the source of considerable technical difficulties in [7J . 
As the reader will see, the proof of Theorem U] is much more straightforward than 
that of Theorem 5 in [7J . 

A direct proof of Proposition [5] is provided in Sectional However, it could also 
be deduced from Perkins' disintegration theorem together with Le Gall's excursion 
representation of Dawson- Watanabe processes. For this, one must use the fact that 
the empirical measure for an exchangeable system of particles fully determines the 
law of the individual trajectories. 

In light of Theorem Q] and the above discussion, one can wonder if, conversely, 
Theorem Q] could not be extended to Beta-coalescents. Interestingly, we have rather 
strong evidence that this cannot be true. Assuming that the construction was valid 
for a Beta-coalescent, we would conclude that the number of blocks lost at every 
collision would be an i.i.d. sequence. Indeed, from results in [10] . one can see that 
the tree spanned by vertices at a deterministic level (the "reduced tree") forms a 
time-change of a Galton- Watson process. However, in that case the distribution of 
the number of blocks involved at every collision must agree with the distribution 
in equation (10) of [6]. Simulations indicate that this is not the case for the final 
collision. However, we are not aware of a rigorous proof that this cannot be true, 
although there are also some theoretical arguments in that direction. We note that 
the essential reason for which the proof of Theorem Q] breaks down in the case a < 2 
is that Brownian motion possesses an extra symmetry property (reflection about 
a horizontal line) compared to the height process of other stable continuous-state 
branching processes. 

2.3 More on duality. 

Theorem [T] brings to mind another result which relates Kingman's coalescent to the 
behaviour of Yule processes. (As the reader will see in section [H Yule processes are 
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indeed embedded in the construction of Theorem [T]) . Consider the jump chain of 
Kingman's coalescent i.e., the chain X(l), . . . , X(n), . . . where X(n) is the element 
of the n th simplex given by the block frequencies of Kingman's coalescent when 
it has n blocks. In [9], Bertoin and Goldschmidt show that this chain X can be 
related to the fragmentation chain obtained by taking a Yule process (Y t ,t > 0) 
and conditioning on {W = w > 0}, where W = ]imt- too e _t Yi almost surely. This 
fragmentation process (Gt,t > 0) is defined by considering how the mass w is 
being split between the children of an individual when it dies. More precisely, 
the two children of this individual each give rise to their own independent Yule 
process , say. To each one we can associate the corresponding random 

variable and W^. Note that if r is the time of the split, we must have 

W = e~ T W^ + e~ T W^ 2 \ Define the fragmentation process (Gt,t > 0) by saying 
that at time r, the mass W splits into e~ T W^' and e~ T W^ 2 \ and so on as time 
evolves. If we do not condition on {W = w}, then G is a fragmentation in the 
sense that fragments evolve independently with individuals fragments behaving as 
a rescaled copy of the original process. This is what is usually called a homogeneous 
fragmentation process (see, e.g., [8j Section 3.1]) but note that here the total mass 
is random. It is then natural to ask what happens when we condition on {W = w} 
to keep the total mass of this fragmentation deterministic. Theorem 3.1 of [9] then 
states that if (Nt,t > 0) is an independent Poisson process, then conditionally on 
{W = w}, 

{G los{1+t) ,t> 0} = {wX(N wt ),t> 0}. (5) 

That is, up to a time-change, the fragmentation G is the time-reversal of the jump 
chain of Kingman's coalescent. It is an open problem to decide whether a similar 
representation exists for the case of Beta-coalescents. We indicate that Christina 
Goldschmidt [14] has recently computed explicitly the Martin boundary of the 
continuous Galton- Watson process (Yt,t > 0) associated with the reduced tree 
at level 1 of a stable continuum random tree with index a. This can be used 
to describe explicitly the behaviour of Y^ w \ where denotes the process Y 

conditioned on {W = w}, (here again, W = lim^oo e- xt Y t and A = l/(a - 1)). 

In particular, from her result it follows that the transition rates of Y^ w ' are not 
independent of the current state of Y ^ . Since it turns out that W = Z\ almost 
surely, {Z\ being the quantity of local time of the stable CRT at level 1), it follows 
that conditionally on {Z\ = z}, the number of children of individuals in the tree 
is not i.i.d. In particular, the above objection against the extension of Theorem [T] 
does not hold here. 
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3 Brownian construction of Kingman's coa- 
lescent. 



3.1 Proof of Theorem [TJ. 

Proof. The proof we give is based on a few simple calculations with excursion 
theory. For basic definitions and facts, we refer the reader to Chapter XII of |22j . 

Fix < t < 1, and recall our notation Z x = L(C, %) for all x > 0, where L(t, x) 
is the joint local time of a Brownian excursion (B s ,s > 0) conditioned to exceed 
level 1. Define a filtration 

{g u = a(Z s ,u < s < l)}o<u<i- (6) 

In parallel, define a family of cr-algebras {7~C u }o<u<i by putting 

H u = o-(B a u {s) ,s> 0) (7) 

where for all s > 0: 

a u (s) := inf jt > : J 1 {Bt>u} oIt > s 

In words 7i u contains all the information about the trajectory above level u. It is 
a tedious but easy exercise to see that {Hi_(}o<ki is a filtration and 7i u 5 Gu for 
all < u < 1. Let 5 > 0, and let u = 1 — £. Define the event Akj(S) as follows. 
Call N the number of excursions above u which reach 1, and let (ei, . . . e^v) be 
these excursions with an order given by their height (as usual, we view each as 
an element of O*). Let 1 < fe, j < N and let 

j4fcj(<5) := {excursions and ej have coalesced at level u — 5}, 

where we say that and ej have coalesced at level u — 5 if and are part of 
a single excursion above level u — 5 which reaches level 1. 

To alleviate notations we also write A(5) when there is no risk of confusion. 
Theorem [T] follows from the following claim: 

F(A kjj (8)\g u ; N > max(i, k)) = ^5 + o(5), a.s. (8) 

Note that to every excursion e, with 1 < i < n, we may associate a starting time 
ti S [0, (] such that the excursion starts at time ti, that is: B(t + ti) = u + ei(t), for 
all t < C( e i)- Assume without loss of generality that tj < tj.- Among all excursions 
below level u, consider the subset of those excursions (e-,i > 1) which separate e k 
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Figure 2: Between levels u and u—du, excursions ej and (which start respectively 
at times tj and tk) coalesce. For this to occur, the excursions below u starting at 
time ti and t' 2 must not reach below u — du. 

and Cj, that is, whose starting time t\ lies strictly in the interval (tj,tk) (see Figure 
[2]). Observe that for e,- and to coalesce by level u — S, it is both necessary and 
sufficient that: 

mfe' i (s)>-6, foralH>l. (9) 

s>0 

Since L(t, u) may only increase at times t such that Bt = u, it follows that 
throughout ej, the local time at level u is constant, let £j be this quantity (thus, 
£j = L(tj,u)). Similarly, let l]~ be the local time at level u throughout excursion 
e^, and let I = \£j — £k\- Then we claim that, conditionally on £j and £k, A(5) has 
probability exp(— £/{28)). Indeed, by excursion theory and elementary properties 
of Brownian motion, we have that, conditionally on 7i, u , 

i>l 

is a Poisson point process on [£j,£k] x ^* with intensity 1^^ ,e k ]}d£ v u (de), where 
d£ is the one-dimensional Lebesgue measure, and u u is Ito's excursion measure 
restricted to negative excursions e G £1* such that inf s> o e(s) > —u. It follows 
that 

i>l 
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is a Poisson point process on [4^j]xM_ with intensity dZlugiLj, ]}®^-{-u<h<o}dh/2h 2 , 
where z := Z u (see, e.g., (2.10) in Chapter XII of [22]). Thus, 

F(A(5)\H u ) = exp(-\£ j -£ k \/(26)). 

Now, observe that, conditionally on Q u , 

i>l 

is a Poisson point process on [0, z] x M + with intensity d£l^< 2 j (8) lif l >i_ u \dh/2h , 
because e$ are precisely the set of excursions above u that reach level 1. Therefore, 
using elementary properties of Poisson processes, conditionally on Q u , and condi- 
tionally on {N > max(j, fc)} (which we must assume if we are to talk about the 
coalescence of ej and e^), then £j and £ k are two independent random variables 
uniformly distributed on [0, z] (this does not depend on the precise intensity of 
the point process; all that is required is that the time-intensity is the Lebesgue 
measure and the space-intensity has no atom). In particular, conditionally on Q u 
and N > max(j, k), £ := \£j — £k\ may be written as: 

e = z\u-v\ 

where U,V are uniform random variables on (0,1) and z = Z u . Putting these 
pieces together, we obtain: 

F(A(S)\g u ; N > max(j, k)) = E[exp(-z|C7 - V\/(25))], a.s. 

1 a a ( z \ x ~y\\ 

axay exp I — — I a.s. 



o Jo 



1 , f v , ( z{y-x) \ 
ay I ax exp — a.s. 

) Jo V 25 ; 

/ zy \ ( zx \ 

dyeXp {-26) J dxeXp {26) ^ 

dyexp f— 7:7) — (exp (-7) — 1) a.s. 



o 







25 J z v r V2o 

= t/ dy[1 - exp ("S) ] a - s - 

= T-^ (1 - exp v- 2 7J ) a ' s - 

= h o(d) a.s. 

z 

as 5 — > 0. This is precisely (|5]l , and so this implies Theorem [TJ □ 
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3.2 Proof of Proposition [5] 

The proof of Proposition [5] involves similar ideas to the proof of Theorem Q] above. 
However, the calculations become, somewhat surprisingly, slightly more compli- 
cated. In particular, there is a remarkable cancellation towards the end, which 
illustrates the following fact. Roughly speaking, we try to compute the rate at 
which the /c th highest excursion at a current level u > splits to give the k th and 
j th highest excursions, with j > k. We are trying to show that this rate is A/Z u 
and is in particular independent of k and j. This may seem hard to believe at first: 
when splitting, it is easier to create an excursion of small size rather than large. 
However, the excursion heights accumulate near zero and therefore creating an 
interval whose size falls exactly between the sizes of the nth smallest and n + 1th 
smallest intervals also becomes harder when n tends to infinity. These two effects 
exactly compensate each other and imply the aforementioned result! 

Fix k,j > 1 and assume that k < j. Recall that here (B s , s < t\) is no longer 
a single Brownian excursion but a collection of excursions which accumulate one 
unit of local time at level 0. We shall still denote by Zt = L(t±, t) the total local 
time accumulated at level t by (B s ,0 < s < r±). Recall that, by the Ray-Knight 
theorem (Theorem (2.3) in Chapter XI of [22]), (Zt,t > 0) is the Feller diffusion, 
defined by: 

Uz t = sfzldWt 

\Zq = 1 

where (Wt,t > 0) is a standard Brownian motion. Let 

Ku =a(B 1 u {s) ,s > 0) (10) 

where for all s > 0: 

7 u (s) := inf jt > : J l { B T <u}dr > s\ . 

In words 1C U contains all the information about the trajectory of B below level u. 
Just as in the proof of Theorem [TJ {/Ct}t>o is a filtration. We wish to compute 
the rate at which the j th excursion "looks down" on the fc th excursion to adopt its 
label. We claim that, independently of k and j, and independently of Tit (defined 
in (|T|)) the infinitesimal rate at which this happens is: 

4 

rate (j looks down on k) = — . (11) 

Zt 

Here and in what follows, the event that j looks down on k means the following 
thing. Define a process {d(t),t > 0}j>i as follows. Initially, {£«(0)}i>i are i.i.d 
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uniform (0, 1) random variables. We think of £i(0) as the label of the i th highest 
excursion of (B s ,0 < s < n) above 0. As time evolves, £i(t) keeps this inter- 
pretation, with £i(t) being the label of the i th highest excursion above level t of 
(B s ,0 < s < ri). The rule of evolution is that when the level t increases, an 
excursion may split into two parts. The highest of these two parts is necessarily 
still ranked i th , while the second part has a larger rank j > i. In that case, we 
say that j looks down on i. When this occurs, £j(t + ) becomes £i(t~), while for 
k > j, what was previously the k th excursion now becomes the (k + l) th excursion: 
thus in that case, we define £k+i(t) = Visually, j adopts the label of i (it 

looks down on it) and all labels corresponding to excursions with a higher rank are 
pushed up by 1. Thus pip tells us that, in the terminology of Donnelly and Kurtz, 
{£i(t)}i>i is a lookdown process with infinitesimal rate 4/Zt at time t, condition- 
ally on (Zt,t > 0). We do not need the full definition of these processes, which is 
somewhat complicated, and can be found, for instance, in Chapter V of [12]. Thus 
(jlip is exactly the Brownian analogue of Theorem 14 in [7] , where further details 
on this approach can be found. In particular, Proposition [5] follows directly from 
pip and Lemma 5.6 in |12| . 

We give a direct proof of pip based on calculations with Ito's excursion mea- 
sure. Fix a small 5 > 0. For all t > 0, let (et,i)i>i be the excursions of B above 
level t ordered by their height (recall that e^j G for all t > 0, i > 1). For k > 1, 
let 

Mfc(i) = supei. fc (s). 

s>0 

A simple analysis of the construction shows that in time 5 (that is, between times 
s and t := s + S), the j th excursion looks down on the /c th excursion if there is a 
local minimum located within the k th highest excursion at level s<m<t = s + 5 
such that of the two branches going out of this local minimum, one has a height 
equal to M)-(t~) = Mk(t) and the other has a height /i new such that 

Mj (r ) < h acw < M j+1 (*-). (12) 

Let (e/j(x),0 < k < Q be this excursion. Given {M^(t) = m}, has the law of 
an Ito excursion conditioned on the event: {sup xg ( ^ &k{ x ) = m }- F° r > 0, let 
Q[m) denote the probability measure on C([0, oo)) defined by the conditioning of 
v given that {sup s>0 e(s) = m}. Define: 

T} = ini{x > : e k (x) > 5}, 

Tm = the unique x > such that ek(x) = Mk(t) = m, 
Tf = inf{x > T M ■ e k (x) < 5}. 
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Figure 3: Events A x<y on the left and B x ^ y on the right 



We note that (fT2|) occurs if and only if the event AM j+1 ,Mj U BM j+1 ,Mj occurs, 
where 

A XjV := I 3Tg < s <T M : x < sup e k (r) < y and e k (s) < 5 > 
I re[Ti,s] J 

and 

B x , y ■= { x < sup e k (r) <y>. 
{ re[Tf,Cl J 

For future reference, we also let A = Am ]+1 ,m 3 and B = BM j+1 ,M r 
We claim that 

F(A\JC t )=F(B\lC t ) = ^- + o(5), a.s. (13) 
At 

and that P(A n B\K t ) = 0(5 2 ), a.s., so that, up to terms of order o(<5), P(A U 
-B|/C() = 45/Zt, a.s., as required in (jlip . By time-reversal symmetry of the Brow- 
nian excursion, it is clear that P(yl|/Ct) = ¥(B\ICt) so we only do the calculations 
in the case B. 

The key idea to compute P(JB) is to use the Markov property at time T?. (In 
fact an additional step is required, since under this conditioning e k is not exactly 
Markov). Let m > and let < x < y < m be two fixed levels (we think of 
y = Mj and x = Mj + \). We claim that P(B) can be computed as follows: 

Q {m) (B x ^ y )= 6 -[l-^+o{5) (14) 
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Indeed, consider the unconditional problem first. Recall that the part of the tra- 
jectory of a Brownian excursion that comes after Tr 1 is simply a Brownian motion 
(started from 5) run until it hits zero. Let r/ > and let M be the event that, 
starting from 5, it reaches level m and from there returns to zero without reaching 
further than m + rj. Starting from 5, the probability of reaching m before hitting 
is 4>s( m ) = S/m by the standard gambler's ruin probability. Once there, the 
probability of returning to before hitting m + r\ is r//(m + rf). Thus 

F 5 (M) = (f) S (m)r)/m + o(rf) (15) 

On the other hand, the event that, started from 5, both M and B x<y occur can 
be expressed, by using the strong Markov property at time Tg , as: 

P 5 (M n B x y ) = foijn)-?—- o( V ) (16) 

m — o x \ y J 

In the above expression, the second term expresses the fact that the process, having 
reached distance m, only has to return to distance 5 before reaching m + r/. The 
third term corresponds to hitting x before returning to zero, while the fourth 
and last term corresponds to the probability that, having reached level x, the 
process e reaches before y. All these calculations involve the same gambler's 
ruin probability already used above. Combining (|15p with (|16p . we obtain as 
claimed, up to terms of leading order 5: 

QW(B x J = -(l--)+o(6). (17) 

Here the term in o(8) depends only on m = Mk(s). Rewriting (I17p . we find: 

Q(m) {Bxy) = S (l-^j +0 (S). (18) 

Taking the expectation: 

™ ^ K«k-M^))+° (4) 

= S[M(1/Mj(s)) - E(1/M,- +1 ( S ))] + o(5). (19) 

Now, observe that by the Markov property of the underlying continuous random 
tree, we have above level s a Poisson point process of excursions with intensity 
Ito's excursion measure n, run for a local time z = Z s . Thus, the number iV >a of 
excursions higher than some fixed level a > is distributed as: 

N >a = Poisson(z/2a) 
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since n(supe > a) = l/2a. Therefore, 

F(M j (s)>a) = F(N >a >j) 

= i _ _ . . . _ e -^ 2a (z/2a) J -V(i " I)' 

from which it follows that Mj has a density equal to 

-z/2a ( z /2nV 

P(M J e(a,a + da)) = e , {Z ffi 

a{j - 1)! 

This implies: 

E(l/Mj) - E(l/M j+ i) = / -e~ z l 2a [(z/2ay/a(j -1)1- (z/2 a y +1 /a(j)!] da 

io a 
1 1'°° 1 

= 4/ -e~ z/2a (z/2a) j \j/a-z/2a 2 ] da 
J* Jo a 

using the change of variable u = z/2a, we obtain: 

E(l/M,-)-E(l/M j+1 ) = If- /°° e -V — B^-^d«+ /°° e"V— du") 

J- \ Jo z z 2u Jq z J 

ry POO o /"OO 

= — / e -V +1 du - / e-Vdu 

7o *(j - 1)! Jo 

2/(7 + 1)! i! 



* V j! (i - I)' 
2 



The cancellation of terms involving j is quite remarkable. Since z = Z t , we have 
thus proved the claim (fT3j) : 

P(J3) = 2/Z t . 

As explained before, this proves the claim pip . 

4 Proof of Corollary [3] 

Let (-B s ,0 < s < () be a Brownian excursion conditioned to reach 1, and let 
(nt,t > 0) be the realization of Kingman's coalescent described in Theorem [TJ 
Note that by Theorem [TJ one can obtain the frequencies of the blocks of Kt by 
simply considering the excursions of B above level U(t) with a mass proportional 
to the quantity of local time that they accumulate at level 1. Thus, to prove 
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Corollary [3l the first step consists of computing the number of excursions above a 
given level u < 1 which carry a given amount of local time at 1. For < u < 1, 
and x > 0, let N(u,x) be the number of excursions above level 1 — u that carry a 
local time at 1 greater than x. 

Lemma 6. We have the almost sure convergence: 

In 

-^N(u, 2ux) — ► e~ x (20) 
Z\ 

for all x > as u — > 0. 

Proof. To see this, observe that by excursion theory, N(u) := -/V(u, 0) is a Poisson 
random variable with mean Zi- u /(2u). Moreover, conditionally on N(u,0) = 
n, the n excursions that reach level 1 are i.i.d. realizations of the Ito measure 
conditioned to exceed level u. Now, it is well-known that such an excursion, 
upon reaching u, behaves afterwards as Brownian motion killed when returning to 
0, so by standard properties of Brownian motion local time (see Proposition 4.6 
chapter VI in [22]) we have that the amount of local time accumulated by such 
an excursion at level u is exponentially distributed with mean 2u. Therefore, by 
Poisson thinning, conditionally on Z\_ u = z, 

N(u, x) = Poisson ( JL e -*/( 2 «)) . (21) 

Considering excursions with local time greater than 2ux, let 

2ue x 

X u = N(u,2ux). 

Note that, by almost sure continuity of Z s near s = 1, it suffices to prove that 
X u — ► 1 almost surely as u — > 0. Next we recall the following standard Chernoff 
bound for a Poisson random variable Y with parameter ^u: for all h > 0, 

F(Y > n(l + 5)) < e~ htl{1+S) E(e hY ) = exp(/i(e h - 1 - h) - hfiS). 

Since e h - 1 - h ~ h 2 /2 for h — > 0, we see that -A = e h - 1 - h - h5 < for 
sufficiently small h. Thus, using (f2T|) . we are led to the estimate: 



P(|X W - 1| > b\Z x _ u = z)< Ce- Xz ' 2u 
for some A > and some C > 0. Taking the expectation in the above, we get 

¥{\X U -l\>8)< E{Ce~ XZl - u/2u ). 
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For k > 1, let u k = 1/k 2 . Note that for < s < 1, Z s dominates stochastically 
an exponential random variable with mean 2s. Indeed, Z s is greater than the 
local time accumulated by B at level s after T\, the hitting time of 1 by B. Since 
(i?Ti+t)i > 0) has the distribution of a Brownian motion started at 1 killed upon 
hitting zero, we may apply again Proposition 4.6 in Chapter VI of [22], the claim 
follows. We deduce that 

00 00 r oo c -\k 2 x c -x/(2-2/k 2 ) 

£E(Ce-*^/»)<c£jr 2 6 _ 2/fc2 dx 

k — 2 k — 2 

" C ^ Afc 2 (2-2//fc 2 ) + l <0 °' 

Thus by the Borel-Cantelli lemma, we get that 

X Uk — > 1 almost surely as k — > oo. (22) 

Now, to obtain almost sure convergence for other values of u, let k be such that 
u k+i < u < Uk, and consider the process t \— > N(t,2ux),t G Note that, 

for a given size y > say, the difference |iV(«,y) — N(uk,y)\ is bounded by the 
total number of excursions that coalesce during the interval {uk+i, Uf~\. To see this 
observe that u i— > N(u,y) evolves by jumps of size 1 (either two excursions that 
had masses smaller than y coalesce to give birth to an excursion of mass at least 
y and the jump is positive, or two excursions of masses larger than y coalesce, in 
which case N(u,y) decreases by 1). Hence for each coalescence event in (uk+i,Uk] 
the process N(u,y) changes by at most 1. 
Thus we have that for all u G [uk+i,Uk]'- 

\N(u,2ux)-N(u k ,2ux)\ < \N(u k+1 ) - N(u k )\. 

Since N(u,x) is monotone in x, we obtain that for every u £ (u k+ i,Uk\, 

\N(u, 2ux) - N{u k ,2u k x)\ < \N{u k+l ) - N{u k )\ + \N{u k ,2ux) - N(u k , 2u k x)\ 

< \N{u k+1 ) - N{u k )\ + \N(u k ,2u k+1 x) - N{u k ,2u k x)\ 

< 2\N(u k+1 ) - N(u k )\ + \N(u k+u 2u k+1 x) - N(u k , 2u k x) 
= A fc , 

say, where we have set 

A k := 2\N{u k+1 ) - N{u k )\ + \N{u k+1 ,2u k+1 x) - N(u k , 2u k x)\. (23) 
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Multiplying by u and letting X' u = uN(u, 2ux) (so that ([20]) is equivalent to 
X' u — > Z\e~ x /2 almost surely as u — ► for all x > 0), we get: 

sup \X' u — X' Uk \= sup \u(N(u,2ux) — N(u k ,2u k x)) 
ue(u k+1 ,u k ] ue(u k+1 ,u k ] 

+ (u - u k )N{u k ,2u k x)\ 

< sup {uA k + \u - u k \N(u k ,2u k x)} 

ue(u k+1 ,u k ] 

<u k A k + lUk+1 ~ Ukl X' Uk . (24) 
u k 

It is plain that the second term on the right-hand side converges almost surely to 0. 
To see that u k A k — > as well, it suffices to observe that by (|22p applied respectively 
with x = and x > 0, we see that N(u k ) ~ C\k 2 while N(u k ,2u k x) ~ C±k 2 for 
some random C\ % Ci > 0. It follows that both terms in the right-hand side of (|23|) 
are o(k 2 ), i.e., u k A k — ► almost surely. Thus the left-hand side of (|24|) converges 
to as well, and this implies X' u — > Z\e~ x /2 almost surely as ii — > 0. This finishes 
the proof of the lemma. □ 

We trivially obtain from Lemmata result first derived by D. Aldous in [1] (see 
his equation (35)). Let K(t,x) be the number of blocks in Kingman's coalescent 
at time t that are greater than x. 

Lemma 7. For every x > 0, we have the almost sure convergence as t — > 0: 

tK(t,tx/2)^e-*. (25) 

Proof. This is a trivial consequence of Theorem [T] and Lemma [U Indeed, since 
every block at time t corresponds to an excursion above level u with u = U(t) by 
Theorem 1, and since the mass of a block is given by the renormalized amount of 
local time it accumulates at level 1, one may write: 

K(t,tx/2) = N(u,y) 

where u = 1 — U(t) ~ tZi/4 as t — > 0, and y = txZ\/2. Therefore, using the almost 
sure convergence result in Lemma [6] and making the necessary cancellations, we 
obtain the desired estimate (1251). □ 



Corollary [3] follows directly from (125j) . since if B(t) denotes the mass of a 
block randomly chosen among the K(t, 0) blocks present at time t (uniformly at 
random), then F(B(t) > tx/2) = E(K(t, tx/2)/K(t, 0)) — ► e~ x by the Lebesgue 
dominated convergence theorem. It follows that: 

t 
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where E is an exponentially distributed random variable with mean 1. Now, 
observe that the distribution of the frequency F(t) of the block containing 1 is 
nothing but a size-biased version of the law of B(t), that is, for every nonnegative 
Borel function /, we have (denoting B = B{t) and F = F(t)): 

E{f(F)} = E{Bf(B)}. 

By considering for instance for every a > 0, P(i ? < a) = E(Bl^ B<a y), and the 
Lebesgue convergence theorem, we conclude that as t — > 0, 

2t 

where E has a size-biased exponential distribution. That is, E has the distribution 
xe~ x dx which is a Gamma(2) distribution. Equivalently, E is the sum of two stan- 
dard exponential random variables E + E' . This concludes the proof of Corollary 

El 



5 Proof of Theorem 4 



The main ingredient for the proof of Theorem [H is the definition of reduced tree 
associated with our Brownian excursion (B s ,0 < s < Q, and results about the 
multifractal spectrum of the Branching measure of Galton Watson trees due to 
Morters and Shieh[18]. Formally, the reduced tree T can be described by saying 
that (B s ,0 < s < Q encodes a continuum random tree T with a metric d and a 
root o as in [I]. Each vertex z S T has a unique geodesic ^ z : [0, 1] — > T that 
connects it to o, such that if d(o,z) = p, then d(o^ z (t)) = pt for all < t < 1. 
We define T by 

T= |J 7 ,([0,1]). 

z€T:d(z,0)=l 

Informally, T is a continuum random tree obtained from T by taking away ( "prun- 
ing") every vertex whose descendence does not reach distance 1 from the root. 
Equivalently, this is the tree which, at level < u < 1, has as many branches 
as there are excursions above level u that reach level 1. Thus if we let |T(n)| 
be the number of branches of T at level < u < 1, then we have by definition 
|T(n)| = N(u), the number of excursions above u reaching level 1. 

The process |T(u)| is a variant of a process already considered by Neveu and 
Pitman in a seminal paper [20]. The key observation is that: 



{|T(l-e-*)|,t>0}^{y t ,t>0}, (26) 
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where {Yt,t > 0} is a (rate 1) Yule process. This is a continuous-time Galton- 
Watson process where individuals die at rate 1 to give birth to exactly two off- 
springs. In fact, an even stronger property holds. Consider the random tree T 
obtained from T by applying the same exponential time-change as in (|26p . That 
is, for all z G T with d(z,o) = 1, define "f' z (t) = j z (l — e - *). Let 

T:= [J 7 ;([0,oo)). (27) 

zeT:d(z,0)=l 

Then T is a Yule tree. We can thus define the boundary <9T of the tree T by 
taking dT to be the set of rays, i.e., the set of all non-backtracking T- valued paths 
(C(t),t > 0) such that £(i) is at distance t from the root. The boundary dT 
is naturally equipped with a measure \x, called the branching measure, which is 
defined as follows. Thanks to Kesten-Stigum theorem: 

e~*|T(t)| -1+ W (28) 

where W > almost surely. In the Yule case, W is well known to be an exponential 
random variable with mean 1. We next define a metric 5 on dT by declaring that 
for £, £' G dT, <5(£, £') = e~* where t is the time of the "most recent common 
ancestor", that is, t = sup{s > : £(s) = C'( s )}- (^ is the so-called genealogical 
metric on dT). A ball B of radius t for 5 consists of all rays which pass through a 
given vertex z G T, at distance t from the root. The subtree containing z is itself 
a Yule tree and hence the Kesten-Stigum theorem applies to it, let be the 

associated Kesten-Stigum random variable as in (|28|) . Then define 

H{B) = e~ l W{z). 

It is easy to see that [i satisfies the assumptions of Caratheodory's extension the- 
orem (see, e.g., [11], p. 444) and hence defines a finite measure on dT, with total 
mass jii(OT) = W(o). It will also be convenient to define the probability measure 
/,«(•) = M(-)MST). 

Then our key claim is that we can identify (dT,/jP) with Evans' metric space 

(9T,fJ) ^ (S,rj), (29) 

in the sense that we can find a continuous one-to-one map $ : S — > dT such that 
if x G 5 and ( = $(x), 

V {B(x,t)} = ^{B((,l-U(t))}, (30) 

where U(t) is the time-change appearing in Theorem [TJ To do this, we first define 
$ on the integers. For i 6 N, simply define $(i) to be the ray associated with 
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excursion gj, where £1,62,... are the excursions of B above level 1, ordered by 
their respective heights. The definition of $ is easily extended to x € S by taking 
a suitable sequence x n £ N with x n — > x (see the proof of Theorem 5 in [7] for 
further details). 

Having made the identification (|30p , it turns out that Theorem U is now an 
easy application of Theorem 1.2 in [T5] and Theorem 1 in [TO]. These results are 
derived for a discrete-time Galton- Watson tree T, and state the following. Assume 
that T is a discrete-time Galton- Watson tree and that the number of offsprings L 
of an individual is such that for some r > 0, E(exp(rL)) = 00 but E(exp(iL)) < 00 
if t < r, and L is unbounded. Then we have, by Theorem 1.2 in [18] (corresponding 
to the thick part of the spectrum), for all < 9 < a/r: 

dim U € dT : lim sup & = fll = a - rfl, (31) 

L rwoo m n n J 

where a = logE(L) and m = E(L). Moreover when 9 = a/r the above set is 
non-empty (see Lemma 3.3 (ii) in |18j). Assuming further that P(L = 1) > (the 
Schroder case) and letting r = — logP(L = l)/a we have, by Theorem 1 in [19], 
(corresponding to the thin part of the spectrum), for all a < 9 < a(l + 1/t): 



dim (e G OT : hmsup zJggMg^O) = *\ = o( £(i + T ) 
[ n-»oo n J 9 



(32) 



Moreover when ^ = a(l + 1/t), the above set is non-empty almost surely. To use 
these results in our case, let T be the discrete-time tree obtained by sampling T at 
discrete times 1, 2, . . ., so T is a discrete-time Galton- Watson process which belongs 
to the Schroder class. The random variable W = W(o) is of course unchanged so it 
is an exponential random variable with mean 1 (note that W is only equal to half 
the local time Z±, corresponding to a "one-sided" Levy approximation of Z\ from 
below: see, e.g., (1.11) and (1.19) in Chapter VI of [22]). The distribution of L is 
not particularly nice to write down but it is unbounded and we may nonetheless 
identify the parameters a, r and r as follows. Note first that since W is exponential 
with parameter 1, r = 1. Moreover, since 

{e~ t Yt,t > 0} is a martingale 

we obtain m = E(Y"i) = E(L) = e, so a = logE(L) = 1. We obtain r by computing 
¥(L = 1) = e _1 since every individual branches at rate 1, and hence r = 1 as 
well. Using (i30j) together with (f3T|) and (f32l) . it is now straightforward to deduce 
TheoremUJ The details are left to the reader and are similar (in fact, much easier) 
than the proof of Lemma 26 in [7]. Note in particular that, since the construction 
of Kingman's coalescent in Theorem 1 holds at a fixed deterministic level 1, we 
avoid the use of Lemma 24 in [7J . 
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